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ABSTRACT 



Hypervelocity flows for velocities in excess of 1.4 km/sec (Mach 5) require 
very high stagnation temperature to avoid liquefaction. The arc heater wind tunnel 
has been designed to provide such flows. The electric-arc driven wind tunnel can 
develop stagnation temperatures up to 13,000°K which will produce hypervelocity 
flows up to 7 km/sec (earth orbital speed). The nature of the flow, however, is such 
that the high temperature source flow may cause severe gradients at the nozzle exit. 
In order to perform aero-thermodynamic tests the characterization of the flow in the 
test section is required. This paper experimentally determines the stream profiles for 
an arc jet wind tunnel conical nozzle directly from calorimetry and pitot probe 
surveys. 
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Amps 
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Mass flux 
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I. INTRODUCTION 



This report is concerned with developing a procedure by which flow field 
characteristics of very high temperature arcjet flows may be produced directly from 
stream surveys of nozzles. 

Hypervelocity flight in the atmosphere has previously been concerned with 
entry into the atmosphere and descent to the ground. This has been expanded 
during the last decade to include vehicles which will use the upper layers of the 
atmosphere to perform hypervelocity aerodynamic maneuvers and vehicles which 
will fly from the ground into orbit using air breathing powered lift. The design and 
testing of an aeroassisted orbital transfer vehicle (AOTV) and the National Aerospace 
Plane (NASP) will require ground based testing in order to obtain a data base for 
this flight regime. Such data will be used to validate 3-D, reacting real-gas, CFD 
codes that will be used in the design of advanced vehicles where aero- 
thermodynamic effects are important. 

The 60-MW Shuttle Interaction Heating Facility (IHF) was designed to meet the 
testing requirements for development of the thermal protection system of the shuttle 
orbiter vehicle. Its design provided for heat transfer simulation and stagnation- 
region flow simulation for conditions which the shuttle orbiter would encounter. If 
this facility is to be used in aerodynamic testing, the arcjet flow characteristics need to 
be known in order to determine if they provide optimum simulation of the actual 
flow conditions. 

In arcjet operations, extremely high temperatures are generated. For example, 
for the test reported here, the stagnation temperatures are of the order of 600CTK 
(10,800°R). Probes must be highly cooled in order to survive the environment. 
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Conventional techniques using thermocouples, for example, to determine local 
stream energy simply will not work because of these high temperatures. Stream 
tube flow-swallowing techniques can be used up to moderate temperatures, 
however, such enthalpy probes require sharp leading edges; and these can not be 
cooled in the arcjet environment. Moreover, such probes with blunt cooled leading 
edges yield an uncertain stream-tube capture; and resolution and accuracy both have 
large uncertainties. To circumvent these problems for arcjet stream surveys, a simple 
method has been developed that uses small diameter, cooled, hemispherical-shaped 
probes to measure local heat transfer and pitot pressure. It is shown in this report 
that these two measurements yield stream profiles to a high degree of resolution. 
The theoretical method derived here has been formulated into a computer code for 
data reduction of stream surveys. This code can be incorporated into a real-time data 
acquisition system and will provide on-line reduced stream profiles during the 
nozzle survey for immediate assessment. 

PROFILE is a computer code developed to use the 60- MW IHF stream survey 
data for a 13-inch conical nozzle. The PROFILE code directly produces normalized 
velocity, density, mass flux and enthalpy profiles. The derivation of PROFILE uses 
the basic conservation equations and is described in Chapter IV. PROFILE uses the 
heat transfer and stagnation pressure data as inputs and outputs normalized values of 
heat transfer rates, stagnation pressure, velocity, density, mass flux and enthalpy 
with respect to centerline values of the flow. Characteristics of flows for other arcjet 
designs and other nozzles can be readily obtained using heat transfer and stagnation 
pressure survey data in conjunction with the PROFILE code developed here. 
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II. INTERACTION HEATING FACILITY DESCRIPTION 



A schematic diagram and photograph of the IHF arcjet facility are shown in 
Figure 1 and 2, respectively. The major components of the facility include a 
constricted- arc air heater, a 150 MW power supply, two interchangeable hypersonic 
Mach number nozzles, a walk-in test chamber, and ancillary subsystems consisting 
of a steam-ejector vacuum system, cooling water system, and automated data 
acquisition system. The high-power capability of the arc heater (rated to 60 MW) 
provides for full scale subsystem tests in high enthalpy streams in both large scale 
stagnation flows and boundary layer flows at high Reynolds numbers. 

The primary component of the facility is a constricted-arc heater 390 cm 
(153.3 in.) long and 8 cm (3.15 in.) in diameter shown in Figure 3. The rated 
airflow rate is 1 .4 kg/s (3.09 lbs/s) and the rated pressure is 12 atm. Cooling water 
is supplied at 7000 kPa (1000 psi), which provides cooling for a maximum heating 
of 12 kW/cm 2 to the constrictor wall. Eight anode and cathode electrodes, each 
capable of conducting 1000 A, provide a maximum arc current capacity of 8000 A. 
Thus, the capability for continuous operation at high current and high enthalpy 
levels is provided. 

Direct current for operating the arc heater is provided by a 150 MW D.C. 
power supply. Six modules make up the power supply, each consisting of a three- 
phase, full-wave, phase controlled silicon rectifier. One module is capable of 
operating up to a level of 5500 V into an arc load at current levels up to 2700 A. 
The six modules can be connected in any compatible series/parallel arrangement by 
means of remotely operated setup switches. The duty cycle for the power supply is 
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30 minutes on, 30 minutes off. A 150 MW rating is achieved for run durations of 
approximately 15 seconds. 

The facility has two interchangeable nozzles: 1) a semielliptic, Mach number 5 
nozzle with an open-jet, flat-plate test section measuring 80 X 80 cm (31.5 X 31.5 
in.) for boundary-layer surface flow tests and 2) a conical, nozzle for testing large- 
scale stagnation flows in a free jet. The conical nozzle is constructed in 5 sections 
with exit diameters varying from 15.24 cm (6 in.) to 105 cm (41.4 in.). The Mach 
number varies with the nozzle area ratio from 3 to 7.5 depending on the exit 
diameter selected. Both nozzles are water cooled by internal water passages, and the 
throat diameters are 6.03 cm (2.37 in.). Both nozzles are designed to flow into an 
evacuated, walk-in test chamber that can be maintained at a pressure of less than 
0.004 atm (3 mmHg), a pressure that varies with flow rate according to the 
pumping curve of the steam ejector system. Tests have shown that supersonic flow 
with 33 % blockage ratio has been attained in the conical nozzle, corresponding to a 
model diameter of about 60 cm (24 in.). 

An on-line automated data acquisition and computing system is used to 
record arcjet performance and test data. The conversion from analog to digital 
signals is made by means of an integrating digital voltmeter, so that the signals can be 
processed by an on-line computer. Output consists of a real-time display of five 
channels that are monitored during a run, printed and reduced data, and a 
continuous time history data plot. [Ref. 1] 
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III. SURVEY DESCRIPTION 



The design of the IHF allows for the control of two variables: 1) constrictor 
pressure and 2) arc current. The test conditions are uniquely determined by these 
two control settings only. For this survey, the test chamber pressure was 
maintained at a low pressure to allow the measurements to be made in a free jet. No 
attempt was made to control the expansion ratio into the free-jet chamber. In general, 
the flow was as an under-expanded free jet. This type of operation causes the free jet 
to spread and to operate at a slightly higher Mach number relative to the nozzle exit 
plane.The survey measurements close to the nozzle exit plane, then are relatively 
uninfluenced by the free-jet expansion process. Most of the discussion pertains to 
the 7.62-cm (3-inch) station to avoid the complexity of the additional expansion. 
Constrictor pressure was varied for several values of current in order to map the 
operating envelope for stagnation flows using the conical nozzle. 

The probes used to obtain data consisted of a hemispherical Gardon type 
calorimeter and a pitot probe shown in Figures 5 and 6. The probes were 1.58 cm 
(0.625 in) in diameter and were mounted 20.32 cm (8.0 in) apart horizontally. The 
probes were mounted to a transverse mechanism that was positioned to pass 
horizontally through the flow. Sweep measurements were taken at 7.62 cm (3.0 
in), 19.5 cm (7.7 in), 34.8 cm (13.7 in) and 52.6 cm (20.7 in) downstream of the 
nozzle exit. 

Figure 7 shows probes, probe supports and their relationship to the nozzle 
during a test using a 15.24 cm (6.0 in.) exit diameter nozzle. Figure 8 shows 
probes entering the flow. Figure 9 shows probes approximately on centerline of 
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nozzle. All photographs were taken using only radiated light from gas behindthe 
normal shocks on probes and support mechanisms. 

Data consisted of heat transfer rate and stagnation pressure measurements taken 
at 0.212 cm (0.083 in.) intervals across the flow. Typical profiles for heat transfer 
rate and stagnation pressure for each current setting are shown in Figure 10. 
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IV. DERIVATION OF DIRECT METHOD FOR PROFILES 



Determining the stream profiles for an arcjet wind tunnel nozzle requires that the 
heat transfer rate and the stagnation pressure measurements behind a normal shock 
make up the only parameters on which the velocity characteristics depend. Briefly, 
this method is based on the fact that stagnation-point heat transfer is proportional to 
the product pV 3 (approximately) while the pitot pressure is proportional to pV 2 . 
When these two quantities are known, a simultaneous solution gives density and 
velocity directly from the measurements. Once the velocity profile has been 
determined the profiles for enthalpy, mass flux and density can be readily found. 

In the following derivation, free stream flow parameters are those found in front 
of the normal shock and will be designated by subscript one while subscript two 
will denote parameters behind the normal shock. The starting point is the 
momentum and continuity equations in the following form: 



Rearranging the momentum equation and substituting into it the continuity equation 
gives: 




( 1 ) 

( 2 ) 




( 3 ) 



Solving for p2 




( 4 ) 
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Noting that M 2 is much less than 1 , that M 2 = <p( Y, M) and recalling that from state 2 
to stagnation state 2 is along an isentrope then: 




( 1 + (I_1)m 2 ! ) 



Y/(Y-1) 



( 5 ) 



Normalizing equation (5) with centerline values results in 




Pi 




( 6 ) 



Substituting equation (6) into a normalized equation (4) gives: 



pi v i 

'(p,v, 2 ) 

CL 1 CL (7) 

At this point the normalized stagnation pressure behind the shock is related to 
the velocity in front of the shock. The next step introduces the closed form of the 
Fay-Riddell [Ref 2] equation derived in Appendix A., which is a correlation for the 
heat transfer rate to a hemisphere; namely, 



q 0fi4nJ R (h, K) (g) 

where h t is the total enthalpy of the flow and h w is the enthalpy of the wall. The 
units are: h, Btu/lbm; pc, atm; and R, ft. Normalizing with centerline values gives: 
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( 9 ) 




Using the following form of the energy equation 



* i + 





( 10 ) 



and equation (7) in equation (9) while neglecting the enthalpy of the wall compared 
to the total enthalpy of the flow yields: 



Rearranging 





(V’ + Zft,) 



( 11 ) 



2 h, \ 



(?) 



CL 



/ pl 


f v i 1 




V 

< v i J 




(^) 1 
V CL J 




( 2fc,\ 

1+-^ 



1 



CL 



( 12 ) 



Since the velocity in front of the shock is supersonic the last term drops out. 
Changing notation to reflect only centerline and free stream gives: 
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4 



( 4 ) 



CL 




V l ^ 

K (Yi) ri , 

CL J 



( 13 ) 



Squaring equation (13) and dividing by equation (7) removes the density term and 
results in an equation which relates free stream velocity, stagnation pressure and heat 
flux. 




This equation equates the normalized velocity of the free stream to the two measured 
quantities behind the normal shock. 

The assumption made to derive this result are simply that the flow is 
axisymmetric and the enthalpy of the wall is small compared with the total enthalpy 
of the free flow. The wall enthalpy is measured directly from cooling water and is 
found to be approximately 0.3 MJ/kg (134 Btu/lbm) or less than 5 % of total 
enthalpy. 

Once the normalized velocity has been found it is used in equation (7) to find the 
value for the normalized density. 



Pt 



’CL 



Pt 



CL 



CL 



(15) 



Normalized mass flux is obtained from the product of normalized density and 
velocity. 



jnJL 

mf. 



CL 




(16) 



The normalized enthalpy is taken directly from the energy equation where 
gravitational effects are neglected and the static enthalpy is assumed to be much less 
than total enthalpy, which results in the following equation. 
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— = f — f 

h CL V CL J (17) 

Equations (14), (15), (16), and (17) were used in the code PROFILE to calculate the 
characteristic profiles of the arcjet high enthalpy flows using survey data. 
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V. RESULTS 



For high-powered, arc-heated, hypersonic facilities the testing for aerodynamic 
or thermodynamic (heat transfer) studies requires a central core flow with properties 
well characterized. Ideally, of course, constant core flow properties are desirable. 
The present study was undertaken to characterize the flow properties of the 60-MW 
arc heated wind tunnel at NAS A- Ames Research Center. Such flow properties are 
essential to the engineer for interpretation of the results for aero-thermodynamic 
tests. 

Arc heaters operate with a large enthalpy gradients near the walls. These 
gradients cause the flow properties in the central core to vary; and are an important 
concern in simulation in these ground-based test facilities. The degree of this 
variation has not been reported up to the present time for the arcjet facility at NASA- 
Ames Research Center. Enthalpy gradients are inherent in all types of arc heaters 
that operate with cooled walls. For example, the temperature ratio of the wall for a 
typical arc heater varies from 2 - 5 % when the enthalpy level is in the typical vehicle 
reentry range of 18.6 - 41.87 MJ/kg (8000-18000 Btu/lb m ) which in terms of 
velocity equates to 20,000-30,000 ft/sec. The constricted arc heater was developed 
to minimize this gradient by the current density and length to cause the enthalpy 
profile to assume a flatter shape than a free-burning or vortex-stabilized arc heater. 
The theory of Stine-Watson [Ref. 3] shows that for a constricted arc heater the 
enthalpy profile assumes a Bessel function form as opposed to a linear profile for a 
free-burning or vortex- stabilized design (Huels arc) illustrated in Figure 11. 

Results from the profile measurements obtained here can be used to answer the 
question as to whether the property variables (enthalpy, velocity, density and 
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energy flux) have the same profile or shape at the nozzle exit as they do for the 
constrictor. Scant measurements exist for the constrictor. Most of the values for the 
constrictor arc heaters are predicted by the Stine-Watson theory. Elaboration of this 
theory is represented by a computer code called ARCFLO II [Ref. 4], This code 
was developed to predict the performance of high-pressure, high-enthalpy 
constricted-arc heaters. The conservation equations (mass, axial momentum, and 
energy) are written for a steady-state and radiatively participating gas mixture with 
electric current conduction. The governing equations employ the boundary-layer 
approximation and, hence, the pressure in the wall normal direction is assumed 
constant. Auxiliary relations, such as the equation of state, global current continuity, 
and global mass continuity, are used to relate the gas density with the pressure, mass 
flow rate, and mixture properties. 

The profiles of density, velocity, enthalpy and energy flux computed by 
ARCFLO II in the constrictor are compared with experimental survey data obtained 
at the nozzle exit in Figure 12. The difference seen in the density profiles, Figure 
12a, can be accounted for by the increase in density near the wall of the constrictor as 
opposed to the decrease in density in an expanding free jet. The velocity profiles, 
Figure 12b, are very similar with the experimental data showing a slightly flatter 
profile. The enthalpy profiles, Figure 12c, shows for ARCFLO II a central core 
region where the enthalpy drops moderately from the centerline to r/R = 0.4, a 
linearly decreasing region to r/R = 0.95, and then a steep drop near the wall where 
the conduction losses dominate. The ratio of centerline enthalpy to mass average 
enthalpy is 1.6. This low ratio value is typical for wall-stabilized arcs. The survey 
data exhibits a much flatter profile than ARCFLO II with a ratio of centerline 
enthalpy to mass average enthalpy of 1.15. 



13 



The product of the above three profiles is shown in Figure 12d and is the 
energy flux of the airstream. Note that this parameter is proportional to the 
convective heat-transfer rate to a body within the stream and that it is a measure of the 
usable stream area. The flat profile for ARCFLO II implies that stagnation point 
heating to a test body would be constant to within 10 % of the value obtained on the 
centerline flow for r/R up to 0.7. The survey data profile shows a more 
pronounced curve which can be attributed to the decrease in density in the free jet. If 
test chamber static pressure is maintained at the nozzle exit static pressure then the 
density profile would become flatter and thereby improve the energy flux profile. 

Figures 13 - 28 are plots of PROFILE code output data at specified survey test 
points. All test point plots are normalized with respect to the corresponding values 
obtained at the 3-inch downstream positions. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The present study of survey data for the 13-inch IHF nozzle showed: 

That there is good agreement between the ARCFLO II code 
prediction for velocity profile and experimental data; 

Enthalpy profile of experimental data is much flatter than ARCFLO D 
theory predicts, which can be explained by the codes inability to 
account for energy transfer toward constrictor edges. 

It is recommended that: 

Future surveys control test chamber conditions in order to determine 
if improvements in the profile characteristics result; 

Sting mounted heat transfer and stagnation pressure probes be 
developed in order to reduce probe mount interference with 
boundary layer of the flow; 

Survey data be organized into accessible files so that characteristic 
maps for individual nozzles and complete nozzles can be utilized by 
researchers and contractors. 
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APPENDIX A 



DERIVATION OF THE CLOSED FORM OF THE FAY-RIDDELL EQUATION 

by WARREN WINOVICH 

1. The Nusselt number is defined by: 



Nu = 



H , 



(w = wall value) 



(Al) 



2. The heating rate is given by: 



q = H(T S -T W ) 



where 



c = 



P w 



H [h, - h w ) 



c p. c 



h,-K 

c . = ^ = 1 



T -T 

‘ t 1 « 



(by definition of Nusselt number; wall temperature used as reference) 

3. Solving for heating rate in terms of Nusselt number. 



(A2) 



Nu k w ( h , - h w ) 
q = — 



c ,. c 



• _ Nu fiT — i ( K 

[KKK ^ X ^ Pw ^ w 
V 1 L J 



c 



where subscript "e" indicates external flow conditions 



(A3a) 



(A3b) 
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? = 



Nu 

VTe 



V5" a/pTk 



(ft, - ft.) 

/V,, c 



(A3c) 



4. For the laminar boundary layer near the stagnation point: 






therefore 



(A4a) 



Q = 



N 



Vr 






(ft,- K) 

Pr_ c 



(A4b) 



5. Evaluation of velocity gradient at the stagnation point for a hemisphere: 



1 P, 



p t + ~2~g Lu l = P, (momentum equation ) 



u, = 



^ 8 c Pi / p 

-11 - — 

Pr Pt 



(A5a) 



P, 



V 

(Newtonian flow over hemisphere) 




(A5b) 



(A5c) 
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Pr =P, for X = 0 



(A5d) 




6. The heating rate equation becomes: 



Q = 



Nu 



Vr< 



■yj P» 



\/~R 



2 g , 



P, 2 ^ 
pT 



(K-K) 



Pr w c 



1. The Fay and Riddell result for the heating transfer coefficient is: 



Nu 



= 0. 163(Pr) 






VRe 

(equation 63 and 45 Fay and Riddell) 

8. Substituting equation (A7) into equation(A6); using 



P, \ 



v p ',y 



= T,Z S and R = 53.5 



gives: 



q = 0.763(Pr)~ ' 2g c (53. 3) Vpw 



Vr 



q\/~R 



( h , h w) (/> r ) 



0 



^rr V 2 g c (53.3) VpTmT </tX 



(A6) 



(A7) 



(A8a) 



(A8b) 
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9. Using: 



p, = 39.65^- 



(A9a) 



P. 



= 184x 10"’ 




(A9b) 



where 

^ ciact 

^ “ n 

“ formula 

Equation (A9b) is appropriate for temperatures of 6000*K (10000'R) where arcjet 
operate and for reentry and reentry stagnation conditions (V entry > 5000 m/sec ~ 
15000 ft/sec) 



PsPs 



39. 65 x 184 x 10 Pj . 



0.428 



1000 



0.428 



T.Z . 



K, 



Ps 



0.428 

= 9.60 x 10”* — K« 



(A9c) 



P» P» _ T w K, 

p s p s " T W Z W -jTZT K. 

Note that K w = 1 .00 for the expression chosen for viscosity. 

10. The final form for the expression for heating rate becomes: 



(A9d) 



g'V r ^~ 

VM 7 *. - h .) 



0.25 



0.763 ^ 1A -s(r.z.) T, 
— ^/2g c (53.3) 3.10x 10 ^7 

(Pr) ' ( T,Z ,)■ 



0.314 




(A 10a) 
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Simplifying for Pr = 0.71 the equation becomes: 



^fp7(h.~ h w ) 



= {constant } 






(T »2 W ) 



#.*372 



qVR 

^/~pT(h t - h w ) 



= {constant } 



(T.z.) 



f.ltl 



(T.2.) 



#.•372 




Noting that Zw = 1.00 

For Pr = 0.71 the constant in brackets above becomes: 



{constant } = JL2£L </2g c (53.3) 3.10 x 10' 3 = 0.02225 

(Pr) 

and the expression for heating rate is: 



\fp^(h, - h w ) 



{0.02225}] K]* 



T t Z 



#.#372 




(A 10b) 



The term within the braces is only a function of the stagnation state. For usual wind 
tunnel tests the stagnation state is within the limits: 

0.01 < p s < 1 atmosphere 130 < h s < 20,000 Btu/lb m 

For this range of pressure and enthalpy the term in the braces is essentially a 
constant: 



.4 (T,Z 



0.0372 



r.zr 




= 1.873 ± 0.042 
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Therefore the hemisphere formula for heating rate is: 



or: 



'/p7(h.- h w ) 



(0.02225}{1 .873} 



0.0417 



« =0.0417 h.) 

where 

q = Btu/ s ft 2 

p t = atm -» p t (measured) 

R = feet 
h = Btu /lb 

m 



(A 10c) 



Symbols 

Cp Specific heat 

h Entahlpy 

k Thermal Conductivity 

Pr Prandtl Number 

R Radius 

T Temperature 

Z Compressibilty Factor 

p Density 

p. Viscosity 



Subscripts 

e External Conditions 

s Stream Conditions 

t Stagnation Conditions 

w Wall Conditions 
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APPENDIX B 



UNCERTAINTY ANALYSIS FOR DERIVATIONS 



Uncertainty analysis of the heat transfer and stagnation pressure probes is based 

on manufacture's claims for uncertainties. In the case of the heat transfer probe the 

manufacturer list's accuracy of readings between 3 & 4 % of the value recorded. 

For the following analysis the value of 5 % is used. The accuracy of the stagnation 

pressure probe has been found by NASA personnel to be on the order of 1/4 % of 

the recorded value. Using these values the uncertainty of the normalized velocity is 

shown below. 

Starting with equation (14): 

V_\ _ f _0 f Pe ) 

V CLj {icLj[ Pt CLJ 

and taking the logarithmic differential gives: 



Integrating 




d V 


. J_ 1 


(dq_ 


d i CL ) 


1 


Up 


«Pcl) 


V CL 


2 


{ 4 


<icL ) 


4 


l p 


P CL J 



y 2 [( e ( e r)^Ci.] 4 [( e r) P ( e r)/ 7 czJ 

where (e r ) is uncertainty. 

In the worst case where centerline values are opposite in sign of the original 
values results in 
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%¥- = \ [(0.05) - (- 0.05)] - \ [(0.0025) - (- 0.0025)] 

V CL 1 4 



d V 



CL 



L [0 JO] - L [0.005] 
2 4 



so that 



d V 

V CL 



± 5 % 



Applying the velocity uncertainty to the other property equations (15), (16) and 
(17) gives: 

= «) 



d h 



% 



CL 



dp_ 

PcL 




± 10 % 



d mf 

m f C L 




± 15 % 
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APPENDIX C 



COMPUTER PROGRAM FOR DETERMINING PROFILES 

FILE* PR0F2 FORTRAN A1 



PROGRAM PROFILE 

THIS PROGRAM IS DESIGNED TO CALCULATE FLOW PROFILES FOR 
NASA AMES IHF ARCJET WIND TUNNELS. VELOCITY, DENSITY, MASS FLUX 
AND ENTHALPY PROFILES NORMALIZED TO CENTERLINE CONDITIONS 
ARE DEVELOPED. 



DIMENSION QDOT(IOOO), QDOTN(IOOO), PTOT(IOOO), VNOR(IOOO), 

8 RHO(IOOO), MFLUX C 1 000 ) , HTOT(IOOO), DIST(IOOO), 

8 QSQR(IOOO), PTOTN(IOOO) 

REAL VNOR,RHO, MFLUX, HTOT , QSQR, PTOTN, QDOTN 

INTEGER S,I, J,K,L, IMAX, IQDOT , IPTOT , VAL , 1 1 , 1 2 

REAL PMAX , QMAX, VMIN 

CHARACTER x36 G 

CHARACTER XI E,Y,N 

CHARACTER X6 A,B,C 

CHARACTER *45 F 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

THIS SECTION COPIES TEST CONDITIONS INTO OUTPUT FILE 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
READC 5, 1 ) F 
FORMAT ( A45) 

WRITE( 7 , 1 )F 
READC5, 2)G, C 
FORMAT C A36 , A6 ) 

WRITE(7,X)G,C 
DO 10 1=1,3 
READC 5, 1 )F 
WRITE(*,1)F 
WRIT EC 7 , 1 ) F 
0 CONTINUE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
THIS NEXT SECTION READS TEST DATA FOR QDOT AND PTOT. 

FIRST! NORMALIZES VALUES OF QDOT AND PTOT 
SECOND. LOCATES MAXIMUM VALUE TO USE AS CENTERLINE 
THIRDi CALCULATES VALUES FOR VELOCITY, DENSITY, MASS FLUX AND 
ENTHALPY PROFILES NORMALIZED TO CENTERLINE 

NOTEi QDOT AND PTOT DATA SHOULD CONTAIN ONLY ONE ZERO AS LAST 

VALUE. THIS IS CRITERIA USED TO TERMINATE READING OF DATA. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
I MAX = 0 

DO 20 1=1,1000 

READC5, 15)QDOT( I), PTOT(I), DIST(I) 

IFCQDOT Cl) . EQ.OIGOTO 21 
1MAX= IMAX+1 

20 CONTINUE 

21 WRITE(X,X) 'NUMBER OF DATA POINTS READ =',IMAX 
ENDFILEC5) 
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C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C 1 NORMALIZATION OF QDOT AND PTOT 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

QMAX = 0 
PMAX = 0 

DO 200 I = 1 , I MAX 

QMAX = MAX( QMAX , QDOT ( I ) ) 

PMAX = MAX ( PMAX, PTOT(I) ) 

200 CONTINUE 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C 2 DETERMINING LOCATION OF MAX VALUES. THESE WILL BE USED TO 

C ALIGN PROBE VALUES TO CENTERLINE. 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DO 201 1=1 , I MAX 
IQDOT=I 

VAL = DIM ( QMAX, QDOT ( I ) ) 

I F ( VAL . EQ . 0 IGOTO 202 

201 CONTINUE 

202 DO 203 1=1 , IMAX 

IPTOT=I 

VAll=DIM!PMAX,PTOT!I)) 

IFIVALl .EQ.O)GOTO 205 

203 CONTINUE 

205 WRITE(x,x)'I AT QMAX =',IQDOT,* I AT PMAX =',IPTOT 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C 3 CALCULATION OF PROFILES 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

Il=IQDOT+179 
I2=IPTOT+179 

WRITE! X, X ) • 1 1 = ' , 1 1 , • 12 = * , 1 2 

K = 0 

DO 300 I = I QDOT , I 1 
K = K + 1 

QDOTIK K )=ABS(QD0T(I )/QMAX) 

QSQR(K)=QDOTN(K)xx2 

300 DIST(K)=DIST(I)-DIST(IQDOT) 

K = 0 

DO 301 I=IPT0T,I2 
K = K + 1 

301 PTOTN!K)=ABS!PTOT! I)/PMAX) 

DO 302 1=1, K 

VN0R(I) = SQRT(SQRT(QSQR(I)/PTOTN(I))) 

RHO(I) = PTOTNC I)/(l/VN0R(I))xx2 
MFLUX(I) = RHOd)XVNOR(I) 

302 HTOT(I) = VN0R(I)XX2 
GOTO 900 



900 WRITE! 7 , X) ' QMAX = ' , QMAX 
WRITE! 6 , X) 'QMAX = ' , QMAX 
WRITE! 7 , X ) ' PMAX =',PMAX 
WRITE! 6 , X) ' PMAX =',PMAX 

WRITE! 7,X) »==r =============================================== ====r 

8 ================ ' 

WRITE! 7, X) * HEATING STAG VELOCITY DENSITY MFLUX 
8ENTHPY DI ST ' 

WRITE! 7 , X ) ' RATE PRESSURE* 

WRITE!7,X)' ======================================================= 

$================• 

WRITE!7,600)!QDOTN!I),PTOTN!I),VNOR!I),RHO!I),MFLUX!I),HTOT!I),DIS 

8T!I),I=1,K,9) 

600 F0RMAT!6!5X!F5.3)),5X,F6.3) 

STOP 

END 
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OPERATOR 
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Figure 2 Photograph of the 60 MW IHF with Mach 7.5 conical nozzle. 
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CONSTRICTED ARC HEATER 
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Figure 4 Diagram of interchangeable nozzles for the 60 MW IHF. 
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Figure 5 Gardon type calorimeter probe. 
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Figure 6 Stagnation pressure probe - with support mount. 
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Figure 7 Probes and support mounts outside 6-inch nozzle free jet flow. 




Figure 8 Probes entering 6-inch nozzle free jet flow during survey. 
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Figure 9 Probes near centerline of 6-inch nozzle free jet flow. 






(a) Heat transfer rate verses radius at 3-inch 
position from nozzle exit 




(b) Stagnation Pressure verses radius at 3-inch 
position from nozzle exit 



Figure 10 Typical profiles of the heat transfer rate 
and stagnation pressure, 13-inch nozzle (A/A* = 30) 
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Figure 11 Comparison of enthalpy profiles for 8-cm 
constricted arc. 



36 



4 





(a) Density Profile 



(b) Velocity Profile 
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(c) Enthalpy Profile 



(d) Energy Flux Profile 



Figure 12 Comparison of ARCFLO II and survey data profiles of 
8-cm constricted arc. 
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PROFILE 1.0 1 ENTHALPY PROFILE 
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VELOCITY PROFILE 
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ENTHALPY PROFILE 
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VELOCITY PROFILE 1W hTTh>w ENTHALPY PROFILE 
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VELOCITY PROFILE 1.2 T ENTHALPY PROFILE 
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VELOCITY PROFILE 
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VELOCITY PROFILE 
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ENTHALPY PROFILE 
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VELOCITY PROFILE 
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ENTHALPY PROFILE 
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VELOCITY PROFILE M \ ENTHALPY PROFILE 
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1.0 4 enthalpy profile 
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ENTHALPY PROFILE 
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VELOCITY PROFILE 
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VELOCITY PROFILE J ENTHALPY PROFILE 





a 

a. 

e 

• 

1 ~ 

a> 

£ 

o 

a 

.& 

"3 

§5 

0 • 

I- 00 
g 

1 s 

~ s 

s 

2 ^ 

O- u 
X! 

y u 

c o 

•r - *-> 

£-c 

S g 

S JjJ 

a 

£ D - 
u oc 
o t> 

c ^ 

8 o. 

s i 

0 < 

u ^ 

a> o 

N ^ 
N 

C S 

J= .2 

CJ 

= ••5 

1 C 

2 c 



u 

t- © 
fS cj 

o u 

fc ~ 
.§0 § 
to u 



52 



VELOCITY PROFILE 1.6 1 ENTHALPY PROFILE 
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